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The frequency-dependent nonlinear dielectric increment of dipolar fluids in nonpolar fluids is often measured 
by using a stationary relaxation method in which two electric fields are used: The static direct current (DC) field 
of high strength causing the dielectric nonlinearity, and the probing alternating current (AC) field of low strength 
and high frequency. When a nonlinear composite is subjected to a sinusoidal electric field, the electric response 
in the composite will, in general, consist of AC fields at frequencies of higher-order harmonics. Based on the 
Frohlich model, we present a theory to investigate nonlinear AC responses of dipolar fluids containing both 
polarizable monomers and dimers. In the case of monomers only, our theory reproduces the known results. We 
obtain the fundamental, second-, and third-order harmonics of the Frohlich field by performing a perturbation 
expansion. The even-order harmonics are induced by the coupling between the AC and DC fields although the 
system under consideration has a cubic nonlinearity only. The harmonics of the Frohlich field can be affected 
by the field frequency, temperature, dispersion strength and the characteristic frequency of the dipolar fluid, 
as well as the dielectric constant of the nonpolar fluid. The results are found to be in agreement with recent 
experimental observations. 

PACS numbers: 78.20.Bh, 77.22.Ej 
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I. INTRODUCTION 

Electric fields of high strength applied to dipolar fluids 
produce a nonlinearity in the dependence of the polarization 
on the field strength. On the other hand, when molecules 
are placed in an alternating electric field, they display re- 
laxation behavior - for small molecules that occurs typically 
at microwave frequencies. As studying the permittivity of 
molecules, or larger entities such as cells, yields detailed 
information about the structure and interactions in the sys- 
tem, frequency-dependent nonlinear dielectric increment of 
dipolar fluids has recently received much attention, see e. g. 
Refs. IIiIqIII^ and references therein. 



Nonlinear dielectric spectroscopy is a useful method for de- 
tecting intermolecular interactions and studying their kinetics 
due to its exceptional sensitivity to a large range of different 
aggregation processes in fluids. This method has been used, 
for example, to investigate the kinetics of cis-lactams cyclic 
dimerization cholesterol aggregation fl^ in nonpolar 

media, and it has been used to measure the nonlinear dielectric 
increment by using a stationary relaxation method Qj^lH^. 
Other recent and exciting developments in dielectric spec- 
troscopy include, e. g., studies of phase transitions in liquid 
crystals 1 8], the dielectric properties of living cells using high- 
Tc SQUIDS at low frequencies |9], and combining it with mi- 
crofluidic devices lHoll to analyze cells, proteins and bacteria 
noninvasively and without chemically modifying the system. 

The nonlinear dielectric effect in dipolar fluids is not a re- 



sult of a single mechanism, but it is caused by several factors: 
(i) the Langevin saturation of the orientation of the dipole 
due to a strong field, (ii) intermolecular interactions (e. g., 
dipole-dipole and hydrogen bonds) which lead to formation 
of aggregates with a compensated dipole moment, e.g., two 
monomers may form a dimer, A + B ^ AB, and (iii) in- 
tramolecular processes may produce new dipoles, e.g., one 
dimer may dissociate to two monomers, AB —^A + B. 

The first Langevin effect decreases the permittivity of the 
fluid, whereas the other two effects increase it. In fact, for ag- 
gregation and dissociation processes [(ii) and (iii) above] there 
is a dynamic equilibrium which may depend on the external 
field. Moreover, in the case of a field, additional dimers have 
to be divided to reach a new equilibrium due to the activation 
energy of the field. 

When a nonlinear composite is subjected to a sinusoidal 
electric field, the electric response in the composite will, in 
general, consist of alternating current (AC ) fields at frequen- 
cies of hi gher-order harmonics JIU fl2l [l3l Il4l llsL fla. UtII . In 
three previous papers |15, 16, ITj, we studied the nonlinear 
AC responses of colloidal suspensions like electrorheological 
fluids etc. 

In this paper, we present a theory to investigate the non- 
linear AC responses of dipolar fluids. We show how the har- 
monics of the Frohlich field can be affected by the field fre- 
quency, temperature, dispersion strength and the characteristic 
frequency of the dipolar fluid, as well as the dielectric constant 
of the nonpolar fluid in question. Thus, by measuring the non- 
linear AC responses of dipolar fluids, it is possible to investi- 
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gate the frequency-dependent nonlinear dielectric increment. 

This paper is organized as follows. In Sec.|n| we derive the 
expression for the nonlinear dielectric increment of a dipo- 
lar fluid consisting of polarizable dimers and monomers, and 
obtain the desired fundamental, second-, and third-order har- 
monics of the Frohlich field. In Sec. Ull] we numerically cal- 
culate the harmonics of the Frohlich field for a simple case in 
which there are no correlations between the orientations of the 
molecules. This is followed by a discussion and conclusion in 
Sec.Hvl 



II. FORMALISM 

Let us consider a system in which polar molecules (dimers 
and monomers) are dissolved in a nonpolar fluid of dielectric 
constant e^. If one applies an external field, a dimer may dis- 
sociate to two monomers. We denote the dipole moment of a 
monomer by /i, that of a dimer by ^' , and assume that /i' ^ /i. 
In addition, we assume that only monomers and dimers exist 
(i. e., there are no multimers present) and that all dimers have 
the same dipole moment. Furthermore, we assume that the 
time chemical reactions need is short enough to be neglected, 
and that both dimers and monomers have an isotropic polariz- 
ability. 



A. Nonlinear dielectric effect 



At a high field intensity the dependence of the dielectric 
displacement D on the field strength becomes nonhnear fisl 
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D = eoE + 47rx£'^E, 



(1) 



where E is the electric field, eo the field-independent part of 
the permittivity, and x is the (nonlinear) susceptibility. To 
determine D, one resorts to Clausius-Mossotti equation for 
three components 
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where is the dielectric constant of the nonpolar host fluid, 
r] = Ni/{Ni + N2) and i = y/^. Ni and N2 are the 
number densities of dimers and monomers, respectively, fee 
the Boltzmann constant, T the absolute temperature, u; angu- 
lar frequency of the applied field, and ti {T2) the relaxation 
time of dimers (monomers). The polarizabilities of dimers 
and monomers are denoted by ai and a2, respectively. 

The desired field-dependent incremental dielectric constant 
eE is El 
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Then, the nonlinear dielectric effect is characterized by the 
ratio Ae/£;2 |18] 
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Next, we consider the orientation polarization -por of a di- 
electric sphere of volume V containing n\ dimers and ?i2 
monomers with a dipole moment /id (see below) embedded 
in a continuum with dielectric constant eoo (see below). The 
sphere is surrounded by an infinite dielectric with the same 
macroscopic properties as the sphere. 

The average component of the dipole moment in the direc- 
tion of the field due to the dipoles in the sphere is given by 



(Mrf • e) = Vvor = 



-u/k-BT 



(4) 



where e denotes the unit vector in the direction of the external 
field, and X stands for the set of positional and orientational 
variables of all monomers and dimers. Here u is the energy re- 
lated to the dipoles in the sphere, and it consists of four parts: 

1 . The energy of the dipoles in the external field. 

2. The electrostatic interaction energy of the dipoles. 

3. The non-electrostatic interaction energy between the 
dipoles. This interaction is responsible for the short- 
range orientational and positional correlations. 

4. The binding energy (denoted as ub) between two 
"monomer dipoles" in the dimer dipoles. 

In Eq. 0, Mrf is given by 
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with /irf = /i(eoo + 2ef,)/3e6, and = ^'(eoo + 2eb)/3ef,. 
The number of monomers and dimers in the system are de- 
noted by ni and n2- Parameter e^o is the dielectric constant at 
frequencies at which the permanent dipoles (i.e., the orienta- 
tional polarization) cannot follow the changes of the field, but 
in which the atomic and the electronic polarization are still 
the same as in the static field. Therefore, e^o is the dielectric 
constant characteristic for the induced polarization. In prac- 
tice, £00 can be expressed using an expression containing an 
intrinsic dispersion 
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where eoo(O) is the dielectric constant at the high-frequency 
limit, and Ae stands for the dielectric dispersion strength with 
a characteristic frequency ojc- 

The external field in this model is equal to the field within 
the spherical cavity filled with a continuum of dielectric con- 
stant eoo, while the cavity is situated in a dielectric with di- 
electric constant eo. This field is called the Frohlich field Ef, 
as given by the Frohlich model jloll 
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In this equation, the higher-order terms have been omitted. 

With the above definitions and using u as defined in Eq. @, 
we now have 



du 
'dE~F 



= -Mrf-e + /3. 



(7) 



Here /3 = dus/dEp < 0, and it is first assumed to be a 
quantity independent of Ep. Then, we obtain 
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Note the subscript indicates the absence of the field. 

To express the macroscopic saturation behavior in micro- 
scopic quantities, higher derivatives of the average dipole mo- 
ment have to be taken into account: 
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Terms of higher than 3rd order are neglected. 
In addition, we have a general relation 



(Md • e) = Vpor = 
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Using Eq. (|3} and the terms in E^ and E, we obtain 
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x[3((M4)„/y)_5((Mj)2/l/)]. 
On the other hand, we may write 
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Hence the nonlinear dielectric increment Ae/i?^ is explic- 
itly expressed. By setting A^i = 0, /i' = and = 1 (no 
dimers in the system), Eq. (|9j can be reduced to Eq. (7.40) of 
Ref. Li&J where a polar fluid in vacuum is investigated. 

It is worth noting that the expression for Ae/E'^ [Eq. (|9}] 
contains no (3, i.e., the binding energy has no effect on the 
nonlinear dielectric increment. In fact, even if /? = (3{Ef), 
the same result is obtained due to the fact that the binding 
energy is not due to hydrogen bonds, rather than from dipoles. 

Let us choose a simple case for numerical calculations. It 
is assumed that there are no correlations between the orienta- 
tions of the molecules. Thus, we have 
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In this regard, the nonhnear dielectric increment is given by 
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where 



(/) = n'^{-5VNi + 5N2n2 - 2N2) 

+2/iV'^(^i"i iV2"2 - iV2 - 5FiViiV2) 

+^i'^{~WNl + SiVini - 2iVi). 

B. Nonlinear AC responses 

Next, we perform a perturbation expansion 1I21I1 to investi- 
gate the nonlinear AC response of the dipolar fluid. In this 
case, the applied field is E = Edc + Eac(ti'), where E^c 
denotes the DC field, and Eac(w) the sinusoidal AC field, 
namely Eac(w) = Eac sin((jjt). As a result, the Frohlich field 
is expressed in the terms of harmonics lfl4l[T5lfT7ll22il as 

Ef = E^p"^ + E^ sin(wt) + 

+E2U, cos(2wt) + E:i^ sin(3wt) + • • • , (11) 

where the DC component {E^p'^^) and the fundamental, 
second-, and third- order harmonics {E^, E2U1 and E^^) are 
analytically given by 
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FIG. 1: The fundamental, second- and third-order harmonics of the 
Frohlich field as a function of the field frequency for three different 
temperatures. 
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In the above derivation, we have used two identities, 

siT?{u}t) = [l-cos(2wt)]/2andsin^(wt) = (3/4) sin(a't) - 
(1/4) sin(3(jji). It is worth noting that the Frohlich field 
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FIG. 2: Same as Fig.Q but for three different £b (dielectric constant 
of the host fluid ). 



[Eq. il in is a superposition of both odd- and even-order har- 
monics, even though there is initially a cubic nonlinearity only 
[see, Eq. Q]. Actually, the occurrence of both odd- and even- 
order harmonics is simply due to the coupling between the DC 
and the AC fields 1I22II . In Eq. il 1> . we have omitted higher- 
order harmonics (e. g., fourth-order, fifth-order, and so on). 
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FIG. 3: Same as Fig.Q but for three different Ae. 



III. NUMERICAL RESULTS 

Based on Eqs. (I10> . and (II 1> . we are now in a position 
to perform numerical calculations. Without loss of generality, 
we choose the following parameters: V = 5.0 x 10^ cm"^, 
711 = 0.25 X 10l^ n2 = 1.0 X 10^^, ^ = 2.5D, fi' = OMD, 
Ti = 0.018 ns, T2 = 0.02 ns, ai ^ 8 x 10"25^.jj^3^ ^ 
2 X 10~25 cm3, T = 300 K, et = 40, eoo(O) = 2.5, Ae = 10, 
ujc = 10^ rad/s, i?ac = 1 V/cm, and E^c = 10^ V/cm, unless 
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FIG. 4: Same as Fig.Q but for three different i 



otherwise mentioned. Here D is Debye unit (10 e.s.u. of 
electric moment). 

FigureQdisplays the fundamental, second- and third-order 
harmonics of the Frohlich field as a function of the field fre- 
quency. It is shown that the harmonics of the Frohlich field 
are sensitive to the frequency of the field. Decreasing the 
frequency causes the second- and third-order (fundamental) 
harmonics to increase (decrease). Thus, it is possible to in- 
vestigate the frequency-dependent nonlinear dielectric incre- 
ment (Ae/E'^) of dipolar fluids by detecting the nonlinear AC 
responses. We also find that decreasing the temperature T 
causes the second- and third-order harmonics to increase due 



to the change in the Frohlich field. This result is in qualita- 
tive agreement with the experimental findings by Hellemans 
etal. In their experiment, they measured the nonlinear di- 
electric relaxation spectra for 10-TPEB dissolved in benzene 
at three different temperatures, 279 K, 288 K, and 298 K. In 
addition, it is shown that the temperature has no effect on the 
fundamental harmonics (note that the three curves in Fig.^a) 
overlap). 

In Fig. |2j we investigate the effect of the dielectric con- 
stant of the nonpolar fluid on the harmonics, and find that 
decreasing £{, leads to increasing second- and third- order har- 
monics, but decreasing fundamental harmonics. This is due to 
the change of the Frohlich field. However, decreasing Ae has 
exactly the opposite effect, see Fig.|3l 

Figure |4] shows the effect of the intrinsic characteristic fre- 
quency ujc- From this figure, we can conclude that the charac- 
teristic frequency for the harmonics is strongly dependent on 
the intrinsic characteristic frequency of the dipolar fluid. 



IV. DISCUSSIONS AND CONCLUSIONS 

In the present paper, we have studied a system contain- 
ing polarizable monomers and dimers. Our theory can be 
extended to deal with multimers as well. Here, we have 
studied the fundamental, second-, and third-order harmonics. 
The extension to higher-order harmonics is interesting and the 
present theory can be extended to include them simply by 
keeping the higher-order terms in Eqs. and (^J. 

Throughout this paper, we have discussed the case of weak 
nonlinearity. Thus, for extracting the harmonics, a perturba- 
tion expansion was performed. In the case of strong nonlin- 
earity, this approach is no longer valid, and we might need to 
resort to a self-consistent method I iH E3l . In addition, the 
dipolar fluid considered here is seen as dispersive. In fact, in 
real systems, the nonpolar fluid is also dispersive. The present 
theory can be extended for that case as well. 

To summarize, we have presented a theory to investigate the 
nonlinear AC responses of dipolar fluids based on the Frohlich 
model. It was shown, that our theory reproduces the known re- 
sults of a polar fluid in vacuum 1 18] and is consistent with the 
experimental observations of Hellemans et al. 1 2] for the non- 
linear dielectric relaxation spectra of 10-TPEB dissolved in 
benzene. For the simple case where there are no correlations 
between the orientations of the molecules, we have obtained 
the harmonics of the Frohlich field analytically by perform- 
ing a perturbation expansion approach. It has been found that 
the harmonics of the Frohlich field are affected by the field 
frequency, temperature, dispersion strength, and the charac- 
teristic frequency of the dipolar fluid, as well as the dielec- 
tric constant of the nonpolar host fluid. The results are found 
to be in agreement with experimental observations. Thus, by 
measuring the nonlinear AC responses of dipolar fluids, it is 
possible to investigate the frequency-dependent nonlinear di- 
electric increment. 
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